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We develop an effective theory for heavy baryons and their excited states. The approach is 
based on the contracted 0(8) symmetry recently shown to emerge from QCD for these states in the 
combined large N c and heavy quark limits. The effective theory is based on perturbations about 
this limit; a power counting scheme is developed in which the small parameter is A 1 / 2 where A ~ 
1/Nc , A/niQ (with A being a typical strong interaction scale). We derive the effective Hamiltonian 
for strong interactions at next-to-leading order. The next-to-leading order effective Hamiltonian 
depends on only two parameters beyond the known masses of the nucleon and heavy meson. We 
also show that the effective operators for certain electroweak transitions can be obtained with no 
unknown parameters at next-to-leading order. 
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(\) ' I. INTRODUCTION 

' Recently, it was pointed out that QCD when restricted to the space of heavy baryons (i.e. states with baryon 
CNj , number one containing a single heavy — b or c — quark) has a contracted U (4) Q or, more generally 0(8) symmetry 
which emerges in the combined large N c (the number of colors) and heavy quark limits H- Of course, in the real 
world neither N c nor itiq (the generic heavy quark mass) is infinite. Nevertheless if they are large enough, QCD will 
possess an approximate symmetry. Approximate symmetries can provide significant phenomenological information 
about a system and have had a long and distinguished history in strong interaction physics. Effective field theory 
(EFT) is a powerful tool that can be used to extract this information. In this article we develop an effective theory to 
fNJ ' describe heavy baryons based on the approximate contracted 0(8) symmetry. This is important since this allows one 
t-H , to make testable predictions about the phenomenology of heavy baryons. This paper serves the role of bridging the 
formal developments of Ref. with phenomenological descriptions of heavy baryons. Our ultimate goal is to study 
the spectroscopy and electroweak decays of heavy baryon states which will be done in Ref. || . 

Before discussing how to implement the EFT approach in the context of heavy baryons it is useful to recall how 
Oh! the EFT program works. Effective field theory is useful when there is a separation of scales between the low energy 
phenomena and higher energy physics. The short-distance physics can greatly influence the long-distance physics. 
| However, when working at low momentum one cannot resolve the details of what is happening at short distances. 
, The idea of EFT is to describe the low-momentum phenomena in a manner which is insensitive to the details of 
short-distance physics. Formally, one would begin with the full theory and integrate out all of the short-distance 
degrees of freedom in a functional integral. This yields an effective action to be used in a functional integral over the 
low momentum (long-distance) degrees of freedom. Clearly, the effective action so obtained cannot be expressed as 
the integral over a local Lagrangian. However, all of the nonlocalities occur at the scale of the short-distance degrees 
of freedom which have been integrated. Thus, to good approximation the effective action can be written in terms 
of a local Lagrangian. The effects of the short-distance physics are now contained in the coefficients in the effective 
Lagrangian. 

Clearly such an approach is, at best, approximate. Fortunately, if the scale separation between the long-distance 
scales and the short-distance scales is large the approximation can be very good. Formally, one can develop a 
systematic power-counting expansion where the parameter is the ratio of the typical long-distance scale to the typical 
short-distance scale. As a practical matter one must truncate the effective Lagrangian at some order. The power 
counting, however, ensures that the effects of the neglected terms are formally down from the leading term by some 
power in the ratio of the scales. This, in turn, ensures that the principal effects of the short-distance physics are 
contained in a few coefficients in the effective Lagrangian. In this manner an insensitivity to the details of the short- 
distance physics is achieved. The low momentum phenomenology depends on a few coefficients summarizing the 
effects of the short-distance physics; the manner by which the short-distance physics produces these coefficients is 
irrelevant. 

In principle, all of the coefficients in the effective Lagrangian can be calculated from the underlying physics. Often, 
however, the underlying physics is either unknown or is calculationally intractable. EFT, nevertheless provides a 
conceptual framework for making predictions about the system. The key insight is that one can treat the parameters 
in the effective Lagrangian as phenomenological inputs fit from experiment. Thus, providing there are more pieces of 
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independent data than there are free parameters at a given order, one can make real, albeit approximate predictions. 
Chiral perturbation theory is a familiar paradigm for the EFT program. 

Let us now turn to the problem of heavy baryons. For simplicity we restrict our attention here to the problem of 
isoscalar heavy baryons — the Af> and A c (generically Aq) and their excited states. The generalization to include 
non-isoscalars such as the Ef, or E c raises no new conceptual issues but complicates the analysis. The approach of 
Ref. H begins with a restriction to a heavy baryon subspace of the full QCD Hilbert space, i.e. we consider only the 
space of states with baryon number one and one heavy quark. The physical picture that emerges from the model- 
dependent analysis of Ref. [|]J|] is that low-lying states of the heavy baryon in the combined large N c and heavy quark 
limits can be represented as harmonic collective motion of the "brown muck" (the technical term for the light degrees 
of freedom) against the heavy quark. An equivalent picture in the combined limit is that of the model of Refs. |J Q 
in which the heavy baryons are treated as harmonic collective motions of an ordinary baryon against a heavy meson. 
The symmetry emerges quite simply from this harmonic oscillator picture: the contracted 0(8) algebra is the algebra 
of all harmonic oscillator creation and annihilation operators and all bilinears made from them. 

The key to the analysis of Ref. [^) which enables this picture to emerge is a consistent power-counting scheme. 
Formally, since we are considering the combined large N c and heavy quark limits, it is useful to consider a single 
parameter which characterizes departures from this limit. Thus, we introduce: 

A = A/m Q , 1/N C , (1.1) 

where A is a typical strong interaction scale and toq is the mass of the heavy quark. The effective expansion parameter 
turns out not to be A, but A 1 / 2 J2|. The low- lying (collective) excitation energies of the heavy baryons are of order 
A 1 / 2 . This should be contrasted with excitation energies for internal excitation of the brown muck and the dissociation 
energy of the heavy baryon (into an ordinary baryon and a heavy meson) which are each of order A {§. Thus, provided 
that A 1 / 2 can be considered to be small, a scale separation exists between the low-lying collective states of the heavy 
baryons from all other excitations. This, in turn, implies that EFT methods can be applied to this problem. 

We note at the outset, however, that in practice A 1 / 2 ~ 3 -1 / 2 is not a particularly small number so that the 
scale separation is not particularly large. This implies that one must work at relatively high order in order to get 
relatively imprecise results. Thus for example, even working at next-to-leading order one only expects relative accuracy 
nominally only of order A. Clearly, the approach will be at best semi-quantitative. Moreover, the expansion will not 
converge for all observables. On phcnomenological grounds it is clear that at best the approach will be useful for 
describing the ground state of the Aq and the doublet of the first excited orbital excitation. The second orbitally 
excited states (which have not been observed in either the A;, or A c sectors) is either unbound or very near the 
threshold and clearly beyond the harmonic limit implicit in the effective field theory. The situation is reminiscent of 
the contracted SU(2Np) symmetries seen for light baryons which emerge in the large N c limit [p~3|— p~7|| . In that case 
the large N c limit predicts an infinite tower of low-lying baryons. However, only the nucleon and delta for Nj = 2 (or 
the octet and decuplet for Nf = 3) are well described in the real world. 

The construction of the form of the effective Hamiltonian or Lagrangian for the system is very straightforward. 
Following the standard rules of effective field theory, once the low-energy degrees of freedom are identified, one may 
simply write down the most general Lagrangian built from them which is consistent with the underlying symmetries 
and includes all terms up to a fixed order in the power counting. Similarly, one could add external sources to the 
theory and impose power counting and thereby determine the from of currents. Of course, in order to compute 
observables one needs more than the form of the effective Lagrangian and effective operators — one needs values 
for the coefficients. These coefficients can either be obtained purely phcnomenologically or by matching with the 
underlying theory. 

In this paper we derive the form of the effective Hamiltonian at next-to-lcading order (NLO) directly from the 
structure of QCD and the power-counting rules . We do this formally by directly making a change of variables to 
the collective degrees of freedom. This has the benefit of relating the coefficients in the effective theory directly to 
ground state expectation values of known QCD operators. In this derivation the structure of the Hilbert space plays 
an essential role. In particular, it is necessary to show that in the combined limit, the Hilbert space can be written 
as a product space of collective excitations of order A 1 / 2 which may be interpreted as excitations of the brown muck 
moving coherently against the heavy quark, and non-collective excitation of order A which correspond to excitations 
of the brown muck itself. In addition, as we will show, the excitations of the two spaces are independent up to 
corrections of higher order in A. 

We also derive results about the effective theory of more phenomenological interest. One of the central achievements 
of this paper is to show that at NLO the dynamics depends on only two free parameters in the effective Hamiltonian 
beyond known masses of the nucleon and heavy meson. This is of phenomenological significance in that the small 
value of the expansion parameter suggests that even for quite crude calculations one should work at least at NLO. Had 
there been a large number of free parameters at this order it would have suggested that the scheme would be difficult 
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to implement in practice as data on heavy baryons becomes available. Another important result of this paper is the 
demonstration that certain operators that are not contained in the Hamiltonian but which play a role in describing 
electroweak properties of heavy baryons are completely determined without the need for phenomenological coefficients 
(up to corrections of higher order in A 1 / 2 ) from the known structure of QCD and the counting rules. This increases 
the predictive power of the effective theory by eliminating free parameters. 

As noted above, an expansion in A 1 / 2 is unlikely to be rapidly convergent. Thus it seems sensible to work at relatively 
high order in the expansion. Here we have stopped at NLO for two reasons: one practical and one theoretical. The 
practical reason is simply that beyond NLO the number of parameters grows to the point where it is likely to exceed 
the number of observables which have any reasonable hope of being measured in the foreseeable future. The theoretical 
reason is that the change of variables technique used to derive the effective theory in terms of known operators in 
QCD is straightforward only up to NLO. Beyond this order, ambiguities in the definition of the heavy quark mass 
enter as do effects from the coupling of the low-lowing states described by the effective theory with higher energy 
excitations. Thus, the treatment becomes more subtle and complex beyond this order. Of course, this does not mean 
that there is no valid effective theory beyond NLO. However considerable effort must be undertaken in going beyond 
this order and this effort is simply not justified in view of the practical difficulties mentioned above. 

This paper is organized as follows. In Sec. |n| we follow Ref. @ and briefly review the derivation of the low energy 



collective variables for heavy baryons from QCD with an emphasis on the power counting rules. In Sec. [H, we show 
that the structure of the Hilbert space as a product space emerges from standard large -/V c analysis. The form of 
the effective Hamiltonian (up to NLO) is derived from QCD in Sec. IV. Effective operators for electroweak matrix 
elements are derived in Sec. [v[ Finally a brief summary is given in Sec 
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II. DYNAMICAL VARIABLES FROM COUNTING RULES AND QCD 

In this section we review the counting rules derived from QCD for heavy baryons in the combined large N c and 
heavy quark limits. We will not reproduce the derivations of Ref. in detail. Rather, we will outline the basic 
strategy of the derivation along with the principal physical and mathematical assumptions. We also collect results 
which will be needed later in this paper. 

The first physical point of significance is to restrict our attention to a limited part of the QCD Hilbert space — the 
part with baryon number one and heavy quark number one. Furthermore we will ultimately restrict our attention 
to the lowest- lying states with these quantum numbers. For simplicity, in this section we will only consider a single 
species of heavy quark at a time as the strong interaction does not couple different heavy flavors. We will generalize 
to two heavy flavors when we consider electroweak matrix elements in the following sections. 

QCD is a field theory and as such has an infinite number of degrees of freedom. We can envision making a 
canonical transformation from the original quark and gluon degrees of freedom to a new set of variables (at the 
quantum mechanical level it becomes a unitary transformation). The goal is to find such a transformation which 
has the property that a certain set of the variables (which we denote the collective variables), generate the low-lying 
collective excited states when acting (repeatedly) on the ground state (the Aq). Working explicitly with the collective 
variables one can calculate directly the properties of the low-lying states. Because we are working in the context of 
canonical transformations it is far more natural to work with the Hamiltonian formulation of the underlying quantum 
field theory (QCD) rather than the Lagrangian formulation. 

Unfortunately, there is no straightforward way to generate such a canonical transformation to collective variables 
exactly. However, it is easy to find a set of variables which behaves this way approximately. That is, using the power- 



counting scheme of Eq. (1.1) one can find collective variables which when acting on the ground state (repeatedly) 
generate the low- lying collective states plus a small component of higher excited states — the amount of higher excited 
state admixed is suppressed by powers of X 1 / 2 . 

One canonically conjugate pair of collective variables is the total momentum of the system P, the generator of 
translations (which is well defined in QCD) and its conjugate variable X, the generator of momentum boosts. This 
pair of collective variables, however, does not induce the internal excitations of interest here. 

One might also wish to introduce the momentum of the heavy quark as a collective variable. The expression for 
Pq appropriate in the heavy quark limit is 

P Q = J d 3 xQ\x){~iD)Q(x), (2.1) 

where Q(x) is the heavy quark field, and D is the three-dimensional covariant derivative. In the heavy quark limit it 
is apparent that the variable conjugate to Pq is Xq defined by 
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Xq = / d 3 xQ\x)xQ(x). 



(2.2) 



The definitions for Pq and Xq in Eqs. (2.1), ( |2.2| ) are clearly valid in the heavy quark limit for the subspace of states 
with a heavy quark number of unity; in this limit the heavy quark behaves as in nonrelativistic quantum mechanics. 
In the limit Pq and Xq correspond to the generators of translations and momentum boosts for the heavy quark. 
Away from this limit the concept of boosting or translating the heavy quark separately from the full system is not 
strictly well defined and the Pq and Xq defined above need not be canonically conjugate. However, even away from 
the formal limit the re will exist some canonically conjugate pair of operators which smoothly goes to the variables 
defined in Eqs. (2.1), (2.2) as the limit is approached. Moreover, by standard he avy qu ark counting considerations 
this canonically conjugate pair of operators will differ from those defined in Eqs. ( |2.l[ ), ( |2.2| ) by an amount of order 
C(A/toq) = 0(A). Here we are interested in the theory up t o nc xt -to- leading order, i.e. up to relative order A 1 / 2 . 
Accordingly, to the order at which we work we can take Eqs. ( |2.l| ), (2.2) as unambiguous definitions. 

Unfortunately, the conjugate pairs (X,P) and (Xq,Pq) are not independent. Clearly [P,Xq] ^ since P is a 
generator of translations and translating the system necessarily translates the position of the heavy quark. Thus, we 
can not directly use (Xq,Pq) as our collective variables to generate low-lying states. However, from these two pairs 
we can generate two independent conjugate pairs. In particular we can construct a linear combination of P and Pq, 
which we denote p and a linear combination of X and Xq, which we denote x, such that x and p are conjugate to 
each other and commute with P and X up to corrections at next-to- next-to-leading order (NNLO). Intuitively p can 
be thought of as the generator of translations of the brown muck relative to the heavy quark and x as its conjugate. 
As we will see below, the variables (x, p) do act as approximate collective variables in the sense outlined above. 

The explicit construction of p and x is detailed in Ref. The key inputs to this construction are Poincare 
invariance, heavy quark effective theory, and a decomposition of the QCD Hamiltonian organized via power-counting. 
The decomposition of the QCD Hamiltonian is most easily accomplished by thinking of the heavy baryon as a bound 
state of a heavy meson and a nuclcon. Standard large iV c QCD analysis allows us to identify the interaction energy 
as scaling as N®, while the nuclcon mass tun scales as N c Jl^Jl^]. Standard heavy quark analysis implies that the 
interaction energy is of order uiq while the mass of the heavy meson, itih, is given by 

m H = m Q + O(m Q ) . (2.3) 
Thus, when constrained to the heavy baryon part of the Hilbert space, the QCD Hamiltonian may be decomposed as: 

H = m H + m N + Hint , (2.4) 
where TLi n t is of order A . Poincare invariance implies that 



while the standard heavy quark treatment implies that 



\X; ; -H] = 



(2.5) 



[X Qi ,H] 



Pq, 
H 



(2.6) 



Using Eqs. (2.4), ( |2.5| ), (2-6) along with Eq.(2.3) and the fact that mjv,m# ~ A, the methods of Ref. j2j allows one 
to deduce that 



P = 



niH 



m H + m N 
m N 



1 



P - P Q ) (1 + 0(A)) , 

V -X n ) I (1 + 0(A)) 



(2.7) 



The expressions for x and p differ from the equivalent expressions in Ref. [^| in that we have exploited Eq. (2.3) to 
express quantities in terms of run which is directly physically accessible rather than uiq which is not. Throughout 
the remainder of this paper we will often eliminate ttiq in favor of m# . 

Thus, we have constructed two independent pairs of conjugate variables: [x,p] and (X,P). These two pairs 
naturally arise in the Hamiltonian. However, as we will see in Sec. in describing electroweak operators these pairs 
are not so natural. Instead it is natural to use the positions and momenta of the heavy quarks (Xq,Pq) and the 
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corresponding independent set (Xi, Pi) which correspond intuitively to the position and momenta of the light degrees 
of freedom. Formally they are defined as follows: 



X e = x + X Q , 

P l = P-P Q . (2.8) 

In order to exploit the collective variables x and p to learn about the low-lying sp ectr u m, it is use ful to evaluate 
multiple commutators of x and p with H. Following Ref. || one can again use Eqs. ( |2.3| ), (2.4), ( |2.5| ), ( piil ) and the 
known A scaling of and m#, i.e. m^,mH ~ A, to obtain 

[x k ,[ Xj ,n]] = -Siif/i (1 + 0(A)) , (2.9) 

where the reduced mass, [A, is given by 

= m N m H 
m N + m H 

More generally, it can be seen by the same methods that multiple commutators of x with TL are suppressed by 
multiple powers of A. In particular, if one takes n commutators of components of x with 7i with n even, then the n th 
commutator scales as A" -1 . For example: 



[x m , [xi, [x k , [xj,H]]]] ~ A 
[x , [x n , [x m , [xi, [x k , [Xj,H]]]]]] ~ A 
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(2.11) 



In contrast, standard large N c counting implies that a displacement of the light degrees of freedom relative to the 
heavy degrees of freedom by a distance of order A shifts the Hamiltonian by an amount of order A : 

e -v-y He *P-y _ n „ ( X °) , (2.12) 



where y is a c- number of order A . Equation (2.12) in turn implies that multiple commutators of p with TC are 
generically of order unity: 

\p k ,\pj,K\] ~ A , 

[Pm, [Ph [Pk, - A° , 

[Po, [Pn, [Pm, [Pi, [Pk, [P],'H}}}}}} ~ A° 

(2.13) 

The power-counting rules from QCD are essentially summarized in Eqs. (2.11) and ( p. 13 ). By themselves, how- 
ever, they are not sufficient to develop a fully systematic power-counting scheme for observables without additional 
phcnomcnological input. In particular, there are two possible scenarios in which the power-counting may be realized 
selfconsistently. Intuitively, the 0(8) or symmetric realization corresponds to the case where the effective potential 
between the heavy quark and the brown muck has a minimum at x = 0, while the 0(4) (g> 0(4) or symmetry broken 
realization corresponds to the situation where the effective potential has a minimum at f ^ 0. The 0(4) eg) 0(4) 
realization is labeled symmetry broken since it breaks rotational symmetry in the A — ■> limit. As shown in Ref. |3] 
the 0(8) (symmetric) realization has the generic property that 

x - A 1 / 4 , p - A~ 1/4 (2.14) 

in the sense that typical matrix elements between low-lying states of polynomial operators which contain x nx and p Up 
scale as A^™ x_ ™ p ^ 4 . In contrast, in the 0(4) <g> 0(4) realization, 

x - A , p ~ A~ 1/4 . (2.15) 
In this paper we will assume that the dynamics are such that the system is in the 0(8) (symmetric) realization. 
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III. STRUCTURE OF THE HILBERT SPACE 



In this section we will show that in the combined large N c and heavy quark limits, the Hilbert space can be written 
as a product space composed of collective and intrinsic excitations, with the collective excitations having energies of 
order A 1 / 2 and intrinsic excitations are of order unity. The physical picture of these two types of excitations is quite 
clear — the collective excitations describe coherent motion of the brown muck against the heavy quark while the 
intrinsic excitations are excitations of the brown muck itself. 

The basic strategy in demonstrating that the Hilbert space is of this form is to exploit the fact that our expansion 
parameter is valid in the pure large N c limit where Witten's Hartree picture is well established [flSf] . The key point 
is that while nature is far closer to the heavy quark limit than the large iV c limit, our expansion is formulated in the 
combined large N c and heavy quark limits for N c A/mQ arbitrary. Thus it applies whether one takes the heavy quark 
limit first, the large N c limit first or takes the limits simultaneously. We will deduce the structure of the Hilbert space 
by taking the large N c limit first and then argue that the result goes over smoothly to the combined limit. Ultimately, 
the reason that we can legitimately take the limits in any order stems from the fact that the dynamics of the system 
are determined by the interaction Hamiltonian, Ti. m t = TL — run — mjv- Note that 7ii n t is non-singular in both the 
large N c and heavy quark limits, and therefore we can smoothly take either limit or the combined limit first. 

Now consider taking the large N c limit prior to taking the heavy quark limit. In this case we can use Witten's 
Hartree picture of baryons. Each quark propagates in the mean field of the remaining quarks. Low energy excitations 
are described by the excitations of a single quark in the background of the remaining quarks. In the present case 
we can divide these excitations into two classes — excitations of the heavy quark and excitations of one of the light 
quarks. Of course, for the light quark excitation it is necessary to anti-symmetrize the wave function as there are 
many light quarks. In contrast there is only one heavy quark. In the large N c limit these two classes of excitations 
are orthogonal to each other so that the wave function may be written as a direct product: 

|^)n c ^oo = \h)» \£), (3.1) 

where h ( £ ) represent the heavy (light) quark part of the Hartree wave function. 

The excitation energies for both types of excitations are of order N®. However, the excitation energies of the heavy 

— 1/2 

quark are of order uiq . This is straightforward to see; the heavy quark sits in a potential well of depth and width 
independent of toq. As one takes the large toq limit (after having taken the large N c limit), the heavy quark wave 
function becomes localized at the bottom of the potential well and thus acts like an harmonic oscillator. One sees 
that the characteristic excitation energies of the heavy and light types of motion scale with A according to 

AE h ~ A c V Q 1/2 ~ A 1 / 2 , 

AE e ~ N°m° Q ~ A° . (3.2) 

As noted above, in this picture the excitations of the heavy quark and the light quarks are independent. The total 
excitation energy is simply the sum of the heavy excitation energy and the light excitation energy. 

If we now consider the combined limit with the large N c and heavy quark limits taken simultaneously the picture 
changes. The Hartree approximation is no longer valid. However, the general counting rules remain — they simply 
have to be reinterpreted. Physically what happens in the combined limit is that instead of the heavy quark oscillating 
in the static field of the brown muck, as in the Hartree picture, they oscillate against each other. However, the 
structure of the Hilbert space is not altered in its essence. The space can still be broken up into a product space of 



independent classes of excitations as in Eq. (3.1), but now the two classes of excitations are not the heavy and light 
degrees of freedom as in the Hartree approximation. Rather, they are the collective excitations (between the brown 
muck and the heavy quark) and the intrinsic excitations (of the brown muck) : 

|V)a-o = |C> ® |I) , (3.3) 

where C ( I ) represent the collective (intrinsic) part. Similarly, the excitation energies behave as in Eq. ([T^) but 
with the more general description of collective and intrinsic classes of excitations replacing the heavy and light classes 
of the Hartree approximation: 

AE C ~ A 1 / 2 , 

A£j ~ A . (3.4) 

We see that the structure of the Hilbert space is quite simple in the A — ► limit. In order to understand the 
structure of the effective theory, it is necessary to understand how corrections to this leading behavior emerge. In 
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particular, we will need to quantify the mixing between the collective and intrinsic excitations due to sub-leading 
terms in A. Again, our strategy is to pass to the case where the large N c limit is taken first and we have the Hartree 
approximation as a well-defined tool. Then use the smoothness of the two limits to pull back to the case of A — > 
keeping N c A/rriQ arbitrary. 

It is useful to consider the way l/N c corrections emerge on top of the Hartree picture as these corrections will 
ultimately become corrections in our A counting. We begin by assuming that at finite, but large N c , a Hartree 
picture is a good first approximation. We also make the standard large iV c assumption that corrections to the 
energy and to energy splittings are generically of relative order 1/N C . Moreover, these leading relative corrections are 
independent of toq as the heavy quark only contributes 1 out of N c quarks. Thus, these relative corrections due to 
effects beyond the Hartree approximation are necessarily of order A and do not contribute at NLO (which is relative 
order A 1 / 2 ). It is straightforward to see that 1/N C corrections to the energies imply that Hamiltonian matrix elements 

— 1/2 

between two different states with different intrinsic states are generically of order N c . This can be shown in two 
ways. First, off-diagonal matrix elements between different intrinsic states of order N c yield N^ 1 corrections to 
the energies and we know corrections of this order exist. The second way is to construct an explicit wave function 
using a transition induced by the gluon exchange. The gluon exchange contains two factors of the strong coupling 
and scales as N~ , there is a combinatoric factor of N c as the gluon can connect to any quark; and finally, there is 

— 1/2 

a normalization factor of N c because there are iV c distinct terms in the excited state wave function as any of the 
quarks may be excited. Now we can pass from the large N c limit back to the combined A — > limit (with N c K/mq 
fixed) by exploiting smoothness. We obtain the following important result: 

(C,I|W|C,I / )~A 1 / 2 , (3.5) 

where the state |C,I) is a product state 

|C,I) = |C)®|I). (3.6) 



IV. CONSTRUCTION OF THE EFFECTIVE HAMILTONIAN 



In Sec. we reviewed the power-counting rules for the collective variables x and p in the heavy baryon sector of the 
QCD Hilbcrt space. As discussed in Sec. |, the general principle of effective field theory implies that the form of the 
effective theory is completely fixed by the symmetries of the theory and the power-counting rules. Rather than simply 
invoking this general strategy, it is instructive to derive the effective field theory at next-to-leading order directly from 
QCD through the change in variables discussed in Sec. [n| We do this because it is useful to verify that the proposed 
scaling rules are, indeed, selfconsistent; the derivation gives significant insights into the underlying structure of the 
theory; and finally, it is satisfying to verify explicitly the validity of the general EFT philosophy in this case. 

In the present context, where we have been implementing canonical transformations to coll ectiv e variables, it is 
simpler to work with an e ffe ctive H amil tonian rath er than an effective Lagrangian. From Eq. (2.5) and the power- 
counting rules of Eqs. (2.9), (2.11), (2.13) and (2.14) it is straightforward to deduce the form of the QCD Hamiltonian 
in the heavy baryon sector written in terms of the collective variables up to order A (NLO). It is given by 



Hqcd 



m N 



P 2 



2(m N + m H ) 



El 
2[i 



s(0) 



-(3) , -(4) 

C j,k,l XtX 3 Xl "+" 2—1 C j,k,l,m X i X j X k x l 



j,k,l 



j,k,l,m 



+ E 

3 

0(A 3 / 2 ) 



E-(2) 
c j,k x i Xk 



(4.1) 



where m#, mjv ~ A -1 , and c^ a ' ~ A . From rotational invariance, the are completely symmetric Cartesian tensor 
operators of rank a. Note that while terms up to x 4 contribute at this order, only terms up to p 2 contribute. The 
difference between these two is due to the essentially nonrelativi stic na ture of the k inema tics. Formally, this difference 
is implicit in the different sca ling relations for x and p in Eqs. (2.9), (2.11), (2.13) and (2.14). 

Since the form of Eq. (4.1) was deduced from commutation relations, one cannot immediately conclude that the 
c(°) operators are c-numbers. Rather, at this stage in the derivation we can only deduce that they are operators which 
commute with all of the collective variables: 



] = [X,c^] = [p 



c (a) ] = [P,c (a) ] = 0. 



(4.2) 
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Thus, the tensors are, in general, operators which depend only on the non-collective variables, i.e. the 
are purely intrinsic operators. Note that to the extent that the coefficients c^ a ' in the Hamiltonian of Eq. (4.1) are 
operators as opposed to c-numbers, the Hamiltonian is not strictly an effective Hamiltonian. Rather it is the full 
QCD Hamiltonian written in terms of our collective degrees of freedom. The non-collective degrees of freedom have 
not been integrated out; rather their effects are contained in the non-collective operators c^ a \ 

Non-collective operators might be expected to influence the dynamics of the collective degrees of freedom even 
though they all commute with the collective degrees of freedom. For example, consider the ground state matrix 
element of (c^Xi) 2 . This can be expressed in terms of off-diagonal matrix elements: 



(O^Xi^XiM = £ (OKc^x^acKacKc^ ) x i )|0> + 



E 

Pnc 



(oKcf^l/^cK/WKcf^oio) 



(4.3) 



where |0) is the ground state of the heavy baryon system and the intermediate states are grouped into purely collective 
excitations \ac) (with excitation energies of order A 1 / 2 ) which excite only the collective degrees of freedom leaving 
the brown muck unexcited, and non-collective states \Pnc) (with excitation energies of order A or greater) which 
involve an intrinsic excitation of the brown muck. If the operators in the Hamiltonian had been purely collective as in 
a true effective Hamiltonian (which by construction has all non-collective degrees of freedom integrated out), then the 
sum would not contain the non-collective states \(3nc)- However, the non-collective operators induce transitions 
to the non-collective states. By acting twice — once to go into the non-collective space, and once to come back the 
collective space — the non-collective operators can affect collective space matrix elements of operators contained in 
the Hamiltonian. 

However, from Eq. ( |3.5| ) we know that terms which connect different intrinsic states are suppressed by order A 1 / 2 . 
We can use this to deduce that the contribution of non-collective intermediate states to matrix elements of composite 
operators such as in Eq. (4.3) will always be suppressed by order A relative to the contributions of the collective 
intermediated states and thus only contribute at NNLO. This in turn justifies their omission at NLO. Let us see how 
this works for the decomposition in Eq. From the scaling rules of Eq. ( 2.14 ) it is straightforward to see that 

the contribution from the collective intermediate states goes as A 1 / 2 . A typical term in the sum over states outside 
the purely collective space is 



(OcOilcf^ICI) = |(0cMC)| 2 (Oi| Cj (1) |I) 



A^A. 



(4.4) 



where |0c) and |0i) represent the ground states of the collective and intrinsic subspace s, re spectively, while |C) an d |I ) 
represent excited states, and the notation |C,I) represents a product state as in Eq. (3.6) . The equality in Eq. (4.4) 
follows since Xi only acts on the collective space while c[ only acts on the intrinsic space and the scaling behavior 
follows from Eqs. ( 2.14 ) and ( |3.5|) so that the collective part of the expectation value scales as A 1 / 2 while the intrinsic 
part scales as A. Thus, we see that the contribution of non-collective states is down by a factor of A compared to the 
contribution of the collective intermediate states. 

The 0(A) suppression of the effects of non-collective intermediate states relative t o th e collective ones is generic. 
Although it was explicitly demonstrated for the case of the composite operator in Eq. (4.3), it should be apparent that 
such a suppression occurs for all products of operators cont ained in the Hamiltonian. When calculating Hamiltonian 
matrix elements, the effects of the operators S"^ in Eq. (4.1) inducing mixing with the non-collective states will always 
contribute at relative order A {i.e. at NNLO) or less. Thus, working at NLO one can simply replace the operators 
by their expectation values in the ground state (or more generally in the ground state band — the set of collective 
states built on the ground state of the intrinsic degrees of freedom). Therefore, the operators can be treated as 
c-numbers in the effective theory. 

When the coefficients c^ are treated as operators, rotational invariance does not constrain the allowable forms 
in the Hamiltonian beyond what is given in Eq. (4.1). For example, the term J^cj 1 x i 1S allowable despite the fact 
that x is an L = 1 operator; c^ 1 ) could also be L = 1 operator and couple to x leading to an L = operator in the 
Hamiltonian. Once the c^ operators are restricted to being c-numbers in the effective theory, rotational invariance 
greatly restricts the allowable forms: c-numbers do not transform under spatial rotations. Accordingly, the only 
surviving c^ coefficients are those which multiply rotational scalars. Thus, for example, the cW terms must vanish. 
Treating the &- a ' coefficients as c-numbers and imposing rotational symmetry yields the following form of the effective 
Hamiltonian up to 0(A): 



n 



2(m N +m H ) 
II 

7~^c.m. 



+ {m H + m N ) + c<°) + (tj + \kx 2 ) + ±ax 4 + G(A 3 / 2 ). 



n 



x- 1 



\l/2 



(4.5) 
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where H. c . m . refers to the center of mass motion of the entire system while TL\^ refers to the piece of the Hamiltonian 
whose leading contribution is of order A™ . The coefficients , K and a are of order A and may be expressed as 
ground-state expectation values of , and , respectively: 

(°) =(0|c(°)|0), 



c 

K 



^£(0|cg|0>> 



24 



i E(°lClO). (4.6) 



A few comments about this derivation are in order here. In essence, we did not have to explicitly integrate out the 
non-collective higher-lying degrees of freedom. By making a change of variables to x, p, X and P, we already have 
identified variables which were decoupled from the intrinsic degrees of freedom to the order at which we are working. 
Imagine one explicitly made this change of variables prior to formulating the problem in terms of a functional integral. 
At low order the action could be written as the sum of a collective action plus an intrinsic action without coupling 
between them, and the only effect of integrating out the intrinsic degrees of freedom would be to supply an overall 
constant multiplying the remaining functional integral over the collective degrees of freedom. By making this change 
of variables we make the effect of integrating out the intrinsic degrees of freedom trivial — at least up to NLO. Only 
if we go beyond this order do effects of the the coupling to the intrinsic degrees of freedom play a role, and the effects 
of integrating out the intrinsic degrees of freedom become nontrivial. 

From a phenomenological prospective, it is sufficient to start with Eq. (4.5). The coefficients c^°\ k and a could 



then be determined approximately from fits to experimental data. Given the fact that the masses of the nucleon and 
the heavy meson are known, the effective theory at NLO depends on three parameters — &°\ k and a. The coefficient 
c(°) is determined from the difference in energy between the nucleon-heavy meson threshold and the lowest energy 
state. However, is simply an overall constant and plays no role in the dynamics; the dynamics of the problem 
depends on two parameters at NLO. It is worth noting that at leading order the dynamics depends only on a single 
parameter, k. Thus, going from leading order to next-to-leading costs only an additional parameter, a. 

V. EFFECTIVE OPERATORS 

To make predictions of electroweak current matrix elements between heavy baryon states, one needs more than 
just the effective Hamiltonian. Effective operators are also required. For general composite operators expressed in 
terms of the quark and glue degrees of freedom of QCD, the effective operators can be derived in a functional integral 
form by integrating out the intrinsic degrees of freedom. The most straightforward way to formulate this is as an 
integral over all QCD variables with 5-functions put in to pick out particular values of the collective variables. In 
practice, however, the preceding description is rather formal as it is not generally possible to perform this functional 
integration. At a phenomenological level, one can simply write a form for the operator consistent with symmetries 
and parameterized by some constants and then fit the constants from experimental data. However, if the operators 
of interest can be expressed entirely in terms of the collective degrees of freedom at the QCD level, the functional 
integral approach becomes unnecessary (for the same reason it was unnecessary in the derivation of the collective 
Hamiltonian). At NLO, one can immediately express the effective operators in terms of the same collective degrees 
of freedom. 

To see how this works, consider some operator which at the QCD level can be expressed entirely in terms of the 
collective operators. We seek to find an analogous operator for use in the effective theory whose matrix elements 
in the effective theory reproduces the matrix elements of the QCD operator in the QCD states (up to correction at 
NNLO). Let us start by considering the matrix elements at the QCD level. From the discussion of the structure of 



the Hilbert space at the QCD level in Sec. [II, it is apparent that a low-energy (i.e. with excitation energy of order 
A 1 / 2 ) heavy baryon state in a A 1 / 2 expansion at NLO can be written as: 

\1>) = |C,0i> + A 1 / 2 Y, a c,i|C',I), (5.1) 

C',I^0i 

where |C,I) is defined as a product sta te a s in Eq. ( |3.6| ), and ac.i are constants of order unity. 

The sum in the second term of Eq. ( |5.l|) is restricted to contributions with I ^ Oi due to the product structure of 
the state. If the sum were not restricted, then there would be contributions of the form 
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|C, Or) + A 1 / 2 £ ao',0, |C, Oi ) = |C) ® |0i) + A 1 / 2 £ a C ', 0l |C> ® |0i ) 
c c 

= ^|C) + A 1 / 2 ^ a c ,, 0l |C')^ ® |0i) = |C, Ox) , (5.2) 

where, |C) = (|C) + A 1 / 2 ac, oi C')) is simply a single state in the purely collective space. Thus, the contribution 
to the sum from Eq. (5.2) is of the form of a single product state. If there are contributions of the form as in Eq. (5.2) 
we can replace C by C and thereby eliminate contributions with I = Oi from the sum. Doing this has the advantage 
of removing any possibility of double counting these contributions. 

Now let us evaluate the matrix element of a purely collective operator, Oc, i.e. an operator constructed entirely of 



the collective operators x, p, X and P. From the structure of the state in Eq. ( |5.l| ), the matrix element between two 
low-energy heavy baryon states (given at NLO) is given by 

(^i|e c |^ 2 ) = ((C^OiJ + A 1 / 2 Y, a Ci'Ji <Ci', I J B c I |C 2 ,0i 2 ) + A 1/2 Yl fl c 2 M 2 |C 2 ',I 2 ) 
\ (Vji^Oij / \ c 2 ',i2^0i 2 

= <Ci|e |C a ) + A Y <^ Ocx'a a <V,i 2 (Ci'|e c |C 2 ') = (Ci|e c |C 2 > (1 + 0(A)) . (5.3) 

Ci',C 2 'Ji#0 Il ,I 2 #0t 2 

Thus, up to NNLO corrections the QCD matrix element is simply the matrix element in the collective part of the 
Hilbert space. On the other hand, the effective Hamiltonian we use at this order with the operators cS a ^ replaced by 
c-numbers is designed to correctly describe the collective part of the Hilbert space at this order. 

The preceding argument shows, that in the case of purely collective operators the operators in the effective theory 
are related to the QCD operators in a trivial way up to corrections at NNLO. A purely collective operator at the 
QCD level is of the form 

qQCD (f QCD 5 pQCD ^QCD j pQCD^ ^ 

where the superscript QCD indicates that all of the operators are at the QCD level. The effective operator up to 
order A corrections is then obtained by keeping the same operator structure of ©§ CD , but replacing x^ CD , pQ CD ; 
^qcd anc j pQCD w jth the collective operators x, p, X and P: 

q^s = e^ cv (x,p,x,P). 

This is essentially a nonrenormalization theorem to this order. The effects of integrating out collective degrees of 
freedom only renormalize the effective operators at NNLO. 

In a similar way, all purely non-collective operators at the QCD level, i.e. all operators which commute with all 
of the collective operators, simply become c-numbers in the effective theory with a value equal to the ground state 
expectation value of the analogous QCD operator. These general principles can be used to deduce the form of several 
operators in the effective theory which will prove useful in the calculation of electroweak observables for heavy baryon 
states. 

Let us first consider the operator which generates Lorentz boosts, K . The unitary operator for a Lorentz boost is 
B$ = e . The boost operator is useful, for example, in calculating form factors which involve matrix elements of 
states with different momenta. One can use the boost operator to generate these states from a standard state in its 
rest frame. From Poincare invariance, the boost operator is given at the field theoretic level by 

K = - (XU + UX} = {m N + m H )X + 0(A), (5.4) 



where the second equality follows from Eq. (2.4). Thus, up to NLO at the QCD level K = (m^ + mn)X and from 
the arguments given above, this can immediately be taken over in the effective theory. 

Among the observables which can be studied using the combined expansion are matrix elements of the electroweak 
current J = QjTQi, where subscripts i and j refer to the heavy quark flavors and T is the Dirac structure of the 
left-hand current, r = 7 M (1 — 7 5 ). We will show here that the leading term in the combined heavy quark and large 
N c expansion of the current J can be expressed in terms of the collective variables x, p, X, P. Moreover, we will 
show that the lowest order corrections come at order A and not at A 1 / 2 , so that they contribute only at NNLO. 



10 



The combined expansion of the current J can be done analogously to the pure heavy quark expansion in the Heavy 
Quark Effective Theory (HQET) (20-25 . The idea is to decompose the total heavy quark field into "large" and 



"small" components or "light" and "heavy" fields with the latter being suppressed by a factor of I/tuq relative to 
the "large" component. It is then possible to eliminate the "small" component in terms of the "large" component by 
means of a power series expansion. Near the combined limit considered here, the pure heavy quark expansion breaks 
down. It has been shown in Refs. 26-2^] that the next-to-leading order terms in the pure heavy quark expansion 



of the matrix elements of the current J between the heavy baryon states contain coefficients proportional to A/tuq 
where A = toa q — itlq (with m\ Q being the total mass of a heavy baryon). In the HQET the constant A is treated as 
being much smaller than rriQ, so that the NLO operators are indeed suppressed. The situation is different, however, 
near the combined limit where A is of order m^r, so that A ~ A -1 and A/toq ~ A . Hence, these corrections cannot 
be neglected near the combined limit. 

The appearance of large corrections near the combined limit is a result of the way the total heavy quark field is 
separated into its "large" and "small" components. The definition of these components is necessarily different from 
the analogous decomposition in HQET. Near the combined limit, as in HQET, the heavy quark is not far off-shell since 
its interaction with the brown muck is of order A while its mass is of order A . However, the space-time dependence 
of the heavy quark is much different near the combined limit. Near the pure heavy quark limit the brown muck 
mass is formally much less than the heavy quark mass. As a result, the space-time dependence of the heavy quark is 
essentially determined by the free particle Lagrangian. In HQET, therefore, it is useful to redefine the heavy quark 
field by removing the phase factor of the solution of the free Dirac equation, namely e~ lnlQV Vfi . Near the combined 
limit the heavy quark remains essentially nonrelativistic. However, the space-time dependence of the heavy quark is 
much different from that of a free particle. The phase factor now should contain a contribution due to the brown 
muck mass, which is formally of the same order as the heavy quark mass near the combined limit. This naturally 
leads to the following redefinition of the heavy quark field: 

h$ (y) = e - 4 ( m Q+ m «) P + Q(y) , 

H Q v \y) = e -*(™Q+ m ») ^y^ p_ Q(y) 5 (5.5) 

where is the total 4-velocity of the heavy baryon, the operators P± = (1± ?/)/2 in the rest frame of the heavy 
baryon project out the upper and lower components of the heavy quark field, and y here represents a space-time 
point. As in HQET, the fields and Hq^ satisfy conditions />/ig' = and ^Hq^ — —Hq . Fields liq and 
Hq^ are defined for each value of a heavy baryon 4-velocity as can be seen from Eq. (5.5). The redefinition in 



Eq. (5.5) has the effect of removing both the heavy quark mass, mg, and the nucleon mass, m^, from the dynamics, 
leaving a dynamical scale of order A , i.e., the scale of the interaction Hamiltonian near the combined limit which is 
Hint = H- (m Q + m N ). 

The heavy quark part of the QCD Lagrangian density in terms of the fields hg^ and Hq has the form, 
C Q = Q(ip - m Q )Q = h$ h^> - htf {iv^D a + 2m Q ) H Q v) + h$ (& x - 2m Q ) ffg 



H ( Q V) & ± ) < + m N h%> - flj) 4*) + fffftf - fcWflgO) , ( 5. 6) 



where D± = D — (v v D v )v is the "transverse part" of the covariant derivative D. The heavy quark Lagrangian in 



Eq. (5^) written in terms of the fields hg and Hq differs from its analog in HQET due to the additional term 



proportional to rnjy. This tojv term would naively suggest that the fields hq and Hq^ are both heavy with masses 

(v) 

tojv and rriQ + tun, respectively, which apparently prevents integrating out the "heavy" field Hq . 
This problem is easily resolved if one considers the full Lagrangian density of the system: 

C = C Q + C q + C Y m = Q(ip - m Q )Q + ^ Qj^P ~ m j)lj + ^ym , (5.7) 

3 

where the sum is over all light quarks, and Cym is the Yang-Mills Lagrangian density. The total Lagrangian density 



in Eq. (5.7) should be re-expressed so as to build the brown muck contribution into the heavy degrees of freedom. This 



is completely analogous to the Hamiltonian treatment in Sec. III. In that case it is done by adding and subtracting 
to the Hamiltonian a quantity which is an overall constant m^r and regrouping terms according to their A counting 
scaling. In the Lagrangian formalism this can be accomplished in the following way using a Lorentz covariant operator 
m N Q ^Q: 
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where 



and 



C = (Cq - m N Q i>Q) + (C q + m N Q i<Q + Cym) = C H + Ct , 
Ch = Cq - m N Q iiQ 



(5.8) 



C e = C„ + m N ( 



-YM 



In the rest frame of a heavy baryon operator, Q tfQ is equal to the heavy quark density operator Q^Q and its spatial 
integral equals to one in the Hilbert space of heavy baryons. Thus, in the Hamiltonian the additional term corresponds 
to removing from the light degrees of freedom an overall constant mjv- 

The additional operator proportional to the nuclcon mass can easily be expressed in terms of the fields hg 



and H, 



(«). 



™ n ,/n — m fh( v ) h( v ) ir( v ) u( v ) j_ frWhW iiWijW 
m N Q - m N \ h Q h Q - H Q H Q + H Q h Q - h Q H Q 



This operator cancels the last term in Eq. (5.6) so that the Lagrangian density Ch has the form: 

C H = htf {iv»D») hff - hff (iv»D v + 2m Q ) H Q v) + h$> (ty ± - 2m Q ) + H^ (%> ± ) , (5.9) 



where the condition = is used. Thus, the "large" component hg ' describes a massless field while the "small" 
component Hq has the mass 2uiq. Using equations of motion the "heavy" field Hq can be expressed in terms of 
the "light" field as follows: 



(v) 



1 



2mQ + iv^D, 



-ift ± h 



(?) 



(5.10) 



The relation between fields h^Q and Hq in Eq. (5.10) is identical in form to that in HQET between the "large" and 
"small" components of the heavy quark field Q. 
We can now expand the field Q — e l( - mQ+mN ' )vf ' y ^ (kg +Hq}), the Lagrangian density in Eq. (5.9), and the current 

J = QjTQi in powers of A. This can be done by eliminating fields Hq (using Eq. (5.10)) and expanding the 
numerator in powers of (iv fl D Pi )/2mQ. Thus, the combined expansion of the heavy quark field including terms up to 
order A is, 



Q(y) 



p -i(m, Q +m N )v' J y ll 



1- 



2m 



Q 



h™(y) + 0(A 2 ). 



Using this expansion the Lagrangian density Ch including terms up to order A has the form: 



C H = (<«"D M ) h%> + 



(v) 



1 



2mq 



-h 



(v) 



0(A 2 ), 



(5.11) 



(5.12) 



where g s is a strong coupling constant and G^ u = [-D M , D v ] is the gluon field strength tensor. In a similar way we can 
expand the current J = QjTQi] again keeping terms up to A we get, 



J(y = 0) = &gniW + ^-h^Tmh^ + J±-h% {-if) Th% + 0(A 2 ) . 



(5.13) 



The expressions in Eqs. ( 5.11 ), ( 5.12 ), ( |5.13|) are identical in form to the analogous quantities in HQET. However, 
there is an important difference: these expressions represent the combined expansion in powers of A and not the 1 /tuq 

expansion of HQET. In other words, by defining fields h^Q and Hq' appropriately we were able to resum implicitly 

all the itin /tuq corrections in HQET. 

The correction terms in the combined expansion of the clcctroweak current, Eq. ( 5.13Q , come at order A and not at 
A 1 / 2 . It will be shown in Ref. || that these corrections contribute to the matrix elements at order A as well, i.e., they 
do not induce any anomalously large constants which violate the combined power counting of these operators. Thus, 
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at NLO in our expansion only the leading order operator, h^Th^, contributes to the matrix elements. In addition, 
it will be shown in Ref. || that in the combined limit at NLO there is only a single form factor which parameterizes 
the matrix elements of the current J in the Aq sector. Accordingly, we can consider one Lorentz component of the 
current J, namely the fj, = component. Thus, to completely determine the matrix elements at NLO in the combined 
expansion, we need to know the form of only one operator, namely Hq hg\ 

Let us first consider the flavor conserving operator hg^ ' [y)h^Q [y) . As we are working near the heavy quark limit it 
is legitimate to expand the state in a fock-space decomposition for the number of heavy quarks and antiquarks. For a 
heavy baryon state standard heavy quark analysis implies that the leading fock component has one heavy quark; the 
first sub- leading fock component with two quarks and one antiquark is suppressed by a factor of A/toq ~ A. Thus at 
NLO in A 1 / 2 we can neglect the sub-leading fock component and treat the state as having a single heavy quark. This 
implies that acting in our heavy baryon Hilbert space, the operator Hq {y)fiQ (jf) is 5 3 (Xq — y). Since this is true 
at the QCD level at NLO it immediately holds in the effective theory at that order. 

In order to describe the flavor changing current we need to enlarge the Hilbert space of our effective theory to 
account for electroweak transitions between heavy baryons with heavy quarks of different flavors. The new Hilbert 



space is the direct product of the Hilbert spaces for baryons with a given flavor of a heavy quark described in Sec. III. 
Any state of the enlarged Hilbert space can be written as n/-dimensional vector, where rif is a number of heavy 
flavors. Each component of this vector describes a state with a given flavor. The effective Hamiltonian can be readily 
generalized to include baryons of various heavy quark flavors. It can be written as rif x nf diagonal matrix. For 
example, for two heavy flavors (e.g. charm and bottom) the effective Hamiltonian including terms up to NLO in A 1 / 2 
has the form: 

H eS =( n ^ H M+0(A 3 / 2 ), (5.14) 



where the diagonal entries Hqj and Ti.Q i are the effective Hamiltonians as in Eq. ( |4.5| ). Analogously, QCD operators 
acting in the enlarged Hilbert space become rif x nf matrices. The flavor changing weak transitions can be described 

by off-diagonal matrix elements of an operator hgV {y)^Q.{y), where the subscripts i and j can take integer values 
from 1 to rif. This operator can be written in terms of collective operators only. It follows from arguments analogous 
to the case of the flavor conserving electroweak operator — h^Q\y)h^q\y) — that the flavor changing operator acting 
on the enlarged Hilbert space has the form (for nf = 2): 



WW%W %k-P + ™- (5 - 15) 



The diagonal elements of the operator in Eq. (5.15) correspond to flavor conserving operators, while the off-diagonal 
elements correspond to flavor changing operators. Their action is to instantaneously change the heavy quark flavor 
at the position where a heavy baryon interacts with the current. Since the operator in Eq. ( 5.1 5| ) is written in terms 
of the collective operators only, the corresponding effective operator has an identical form. 
Finally, let us consider the operator defined by 

dt = N l _ 1 Yl J d 3 xql(x)xqi(x), (5.16) 

i— light 

so that di represents the contribution to the baryon dipole moment coming from light quarks. In calculating electro- 
magentic transitions between heavy baryons it is useful to divide the electromagnetic current into pieces coming from 
the heavy and the light quarks and to further divide the light quark contribution into isovector and isoscalar pieces. 
The isoscalar piece is simply the light-quark contribution to the baryon current. Thus, matrix elements of de must be 
calculated to describe electric dipole transitions of heavy baryons. 

At first sight, it might appear to be impossible to express the operator di in terms of collective operators. Although 
di — the collective position operator for light quarks — is intuitively something like the position of the brown muck 
and clearly transforms in the same way as Xi, at the operator level di is manifestly not identical to Xe. Recall that 
by construction m^Xe acts as a generator of boosts for the light degrees of freedom (up to order A corrections). On 
the other hand, the operator dg acts only on the quarks and not on the gluons. Thus m^di cannot boost the gluons 
and therefore does not act as a boost for the light degrees of freedom (which include both quarks and gluons). Thus, 
the operator di is clearly distinct from Xi and there is apparently no simple way to express di entirely in terms of 
collective operators. 
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Remarkably, however, one can show that 

ty'|d<|V> = W\XtW) (5-17) 



for all \ip') in the collective subspace defined in Sec. [II. Thus, when one restricts consideration to the collective 



subspace which dominates the low energy physics, dg is equivalent to Xt. 

To see how this comes about, we start by considering the commutation relations of dg with Pg, Xq and Pq. To 
determine the commutator with P# we exploit the fact that Pg is a generator of translations for the light degrees of 
freedom so that exp(iPg ■ A) q(x) exp(— iPi ■ A) = q(x — A), where q is a generic light quark field. We can use this 
to deduce that 

exp(iPg • A) (d e ^J exp(-iPg ■ A) = exp(iPg ■ A) I — ^— ^ £ / d 3 x qj (x) x q^x) J exp(-iPg ■ A) 

Y c i=light J 

— r ^ I d 3 xql(x-A)x qi (x-K)] = E / d 3 x'qt(x')(x' + K)q i (x') 

i— light / \ i— light 

^— J d 3 x'ql(x')x'q i (x') + A^- I £ j d 3 x' qj {Sf ) q, (f ) ) = dg + A . (5.18) 

z'=light i— light / 



The first equality in Eq. (5.18) follows from the definition of dg, the second equality follows because Pg translates the 
light quark degrees of freedom, the third equality follows from a change of variables, and the final equality follows 
from the definition of d and the fact that 



]T J d 3 x'q\{x') qi {x>) 



i— light 



counts the net number of light quarks, which is iV c — 1 in a heavy baryon containing one heavy quark. Equation (5. IS ) 
implies that 

[dej,Ptk] = tijk , (5.19) 
while the fact that from its definition dg has no heavy quark operators implies that 

[dtj,X Qj ] =0, [d ij ,P Qk }=0. (5.20) 

Thus, the commutation relations of dg with Pg, Xq and Pq are identical to the commutation relations of Xg with Pg, 
Xq and Pq. In Eqs. ( 5.18| ), (5.20) indices i,j = 1,2,3 indicate Cartesian components of Xg, Pg, Xq and Pq. 



Next let us express dg in terms of the various collective and non-collective operators in the problem. Generically, 
one can represent an arbitrary operator O as a power series in our collective operators: 

o = E& n p " Pu n < Xu n n x ^ XQ - ^ 

N i=l,2,3 i=l,2,3 i=l,2,3 i=l,2,3 

with N = (N Pe j , N Pe 2 , N Pe 3 , N Xe 1 , N Xe 2 , N Xe 3 , N Pq i , N Pq a , 7V Pq 3 , 7V Xq x , N Xq 2 , 3 ) specifying a particular term 
in the expansion, and i = 1,2,3 indicating the three Cartesian directions. The coefficients in the expansion, b®, are 
generally non-collective operators which commute with Xg, Pi, Xq and Pq. Expanding one component of dg in this 
form gives 



*i=*i+ E S S> X . II P ^ (5-22) 

Wph.Wp^,^, 4=1,2,3 

where the fact that the (if commutes with Xq or Pq impose the restriction that no nonzero powers of Xq or Pq 



can contribute to the series, while the commutation relation of Eq. (5.19) implies that the only contributing term 



containing any powers of Xg is the linear term and that the coefficient of this term is unity. 
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The known parity and time reversal properties of dg and Pg fix the parity and time reversal properties of the 
non-collective coefficient operators, 5^ Np Np . Since dg is time reversal even and parity odd while Pg is time 
reversal odd and parity even, it follows that: 



Tib 



N P , 



P lb 



,) = (-1)^ 



-JVf 



-Nr 



(-1) 



Nr 



(5.23) 



The significant point about Eq. (5.23) is that every non-collective coefficient operator is odd under PT. From the 
analysis of Sec. Ill, generic states can be written as superpositions of outer product states containing collective and 
intrinsic parts. The collective subspace of the theory is the subspace where the intrinsic wave function is in its ground 
state. Thus, an arbitrary state in the collective subspace can be written as — \ipc) ® |0j) where 0/ indicat es the 
ground state of the intrinsic system. Now consider the matrix element of a typical term in the series of Eq. ( |5~2l ) 
between two collective states: 



M#i,.^.,*P,. II Pt? Pu m = Wc\ II Pe i PU lV>c}<0/|S& t Np Np |0j> , 



Np. 



(5.24) 



i=l,2,3 



i=l,2,3 



wher e we have separated the matrix element into its collective and intrinsic parts. Clearly, the matrix element in 
Eq. (5.24) vanishes: the intrinsic part of the of the matrix element is a diagonal matrix element of a PT odd operator 
in a state of good PT. Thus, within the collective subspace all terms in the series of Eq. ( 5.22] ) vanish except the first 
term . Th is implies that within this subspace all matrix elements of dg are identical to matrix elements of Xg and 
Eq. ( 5.17 ) is established. 

As written, the equivalence expressed in Eq. (5.17) only holds for states in the collective subspace, i.e., states with 
the intrinsic subspace in its ground state. As discussed in Sec. IV, the non-collective operators in the Hamiltonian 
induce components with excited intrinsic states in the physical low-energy states. This in turn spoils the formal 



equivalence in Eq. fl5.17| ) for the physical states. Fortunately, however, it is straightforward to see that in the physical 
states, the matrix elements of dg differ from those of Xg by an amount of order A (or less) and hence may be neglected 
at next-to-leading order in A 1 / 2 . In the first place using reasoning analogous to that leading to Eq. it is easy 



to see that the matrix elements of the b 



Np ei ,Np e2 ,N Pe 



between the intrinsic ground state and its excited states is 



of order A 1 / 2 . On the other hand, the analysis of Sec. fVl implies that admixtures of excited intrinsic states in the 
physical states are suppressed by at least order A 1 / 2 . Combining these two facts one sees that the matrix elements of 
de and Xg in the physical states are equivalent at NLO. 

The equivalence at NLO between the matrix elements of dg and Xg may seem paradoxical. After all, as discussed 
above, the two operators are clearly quite distinct at the QCD level. It would seem that the fact that the operators are 
different would imply that the matrix elements of the two operators should be different in general. The resolution of 
the paradox is quite simple: the matrix elements of the two operators are different in general. They only coincide (up 
to NNLO corrections) for a very special class of states — the collective states with excitation energies of order A 1 / 2 . 
For states with excitation energies of order unity, the matrix elements of the two operators can be quite different. 
Fortunately, we will be using our effective theory to study these low- lying excitations only. Thus, for the purpose of 
the effective theory at NLO, it is legitimate to replace de with Xg. 



VI. SUMMARY 



In this paper we have derived an effective theory which will be used in Ref. (3| to study the spectroscopy and 
electroweak decays of the low-energy states of heavy baryons. The analysis is based largely on the formulation of 
the problem in Ref. P|. A power-counting scheme is developed in A 1 / 2 , where A ~ l/N c ,A/mQ. The kinematical 
variables of the problem A, P, x and p are directly related to QCD operators which are well defined up to corrections 
at relative order A. This formalism is based on standard large N c and heavy quark analysis. 

At the theoretical level, the analysis used in our derivation of the effective theory goes well beyond the analysis of 
Ref. (|] which was largely kinematical and group-theoretical in nature. In our derivation, the structure of the Hilbert 
space for heavy baryons plays a critical role. In particular, as the A — > limit is approached, the heavy baryon 
states can be written as products of collective states (with excitation energies of order A 1 / 2 ) and intrinsic states (with 
excitation energies of order A ). To deduce this form of the Hilbert space we exploited the fact that our formulation 
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as an expansion in A does not depend formally on whether the large N c limit is taken first, the heavy quark limit is 
taken first or they are taken simultaneously. By taking the large N c limit first we can use standard large N c reasoning 
to deduce the structure of the Hilbert space and then map back to the case where the limits are taken concurrently. 
One essential point in this analysis is that matrix elements of operators which couple between different intrinsic states 
are characteristically down by a factor of A 1 / 2 compared to matrix elements of operators with the same intrinsic state. 
Thus, the effect of coupling to an excited intrinsic state and then back to the ground state is typically of order A and 
hence can be neglected at NLO in the A 1 / 2 expansion. 

From the kinematics of the problem, the underlying heavy quark theory, Poincare invariance, and the structure of 
the Hilbert space, we showed that the effective Hamiltonian up to order A can be written in the form of Eq. fl4.5| ). 
At order A 1 / 2 one sees that the excitation spectrum is fixed entirely by a single parameter k, and at order A only one 
additional parameter — a — enters. This is encouraging from a phenomenological perspective since we do not have 
a large number of coefficients to work with at NLO. 

We also derived the certain effective operators at NLO. It turns out that the operators discussed here are precisely 
those needed to compute interesting electroweak transitions of isoscalar heavy baryons at NLO. The key point in the 
derivation is that all of these operators can be written in terms of the collective operators (up to possible NNLO 
corrections) . From this it was easy to show that the effect of the admixture of components only comes in at NNLO so 
that at NLO one can simply replace the QCD level collective operator by an effective operator of the same structure. 
One operator, di, which is effectively the dipole moment of the light quark contribution to the baryon density, was 
somewhat more subtle. Although, as an operator it cannot be written directly in terms of the collective variables, we 
showed that its matrix elements between low-lying states are identical to those of the collective operator Xg. 

The effective theory derived here is suitable for the study of isoscalar heavy baryons at NLO in our A 1 / 2 expansion. 
In Ref. H we will use the effective Hamiltonian of Sec. [V to study the spectroscopy and electroweak decays of A c 
and Ab baryons and their excited states. In order to calculate the semi-leptonic matrix elements we will use relations 
between form- factors analogous to those obtained in HQET |2(i|-|25[ . These matrix elements can then be determined 
using the effective operators discussed in Sec. 0. We will also calculate the radiative decay rates of the first excited 
states of Aq baryons. 
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